TECHNION 96-15 
MZ-TH/96-17 

Spectator Effects in the Heavy Quark Effective Theory 

B. Blok, J. G. Korner, D. Pirjol and J. C. Rojas* 

Department of Physics, Technion- Israel Institute of Technology, Haifa 32000, Israel and 

Johannes- Gutenberg Universitdt, Institut fiir Physik (ThEP), 
Staudinger Weg 7, 55099 Mainz, Germany 

Abstract 

We present a complete analysis of the Heavy Quark Effective Theory Lagrangian at order 
in the leading logarithmic approximation, including effects induced by spectator quarks. At this 
order new correction terms appear in the effective Lagrangian, as four-quark operators containing 
both heavy and light quark fields. We compute the coefficients of these operators in the leading- 
logarithmic approximation, and the matching conditions at one-loop order. Two of these operators 
break the heavy quark spin symmetry, and we estimate their contribution to the hyperfine splitting 
of the heavy mesons. We find that they make a positive contribution to the hyperfine splitting of 
about 10% of the measured splitting in the charm case and of up to 5% in the bottom case. 
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I. INTRODUCTION 



The heavy mass expansion estabhshed itself as a valuable tool for the study of hadrons 
containing one heavy quark 3;|2,y]. This expansion is formulated most naturally in terms 
of an effective theory, the heavy quark effective theory (HQET), which is an approximation 
to QCD with one heavy quark. The corrections to this approximation are controlled by the 
small parameter Ag^/j/m, where Aqcd — 300 MeV is a typical scale of low-energy QCD 
and m is the heavy quark mass. 

The applications of the heavy quark effective theory to physical problems presented so 



ar include correction terms of order J, Sjisl and in some recent instances even 1/m 



.3 



We investigate in this paper the effects induced by spectator quarks in the HQET. 
They appear first at order 1 / and are associated in an effective theory language with 4- 
quark operators containing both heavy- and light-quark fields. These operators mix under 
renormalization with the other operators of dimension six present in the HQET Lagrangian 
already at tree level, and their inclusion changes the form of the renormalized Lagrangian. 

We define the operator basis in Section II and compute their mixing under renormaliza- 
tion. The coefficients of these operators can be obtained from a one-loop matching calcu- 
lation, which is presented in Section HI, together with a leading-log renormalizat ion-group 
running. As an application we estimate in Section IV their total contribution to the hyperfine 
splitting of heavy-light mesons in the framework of the factorization approximation. 

II. DIMENSION-6 OPERATORS 

The Lagrangian of the heavy quark effective theory is written as an expansion in the 
inverse heavy quark mass l|, Q, Q 

= + + + (1) 

The heavy quark field h annihilates static heavy quarks moving with a fixed velocity v and 
satisfies the condition ^{1 + ■^)h = h. The finite-mass effects appear as interaction terms 
in the Lagrangian suppressed by powers of 1/m. The first two correction terms have the 
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well-known expressions j9 



10 



11 



12 



13| 



£i = Ck{^l)h{^Dfh + cM^(^^.vF^"h + U^l)h{lv■Dfh 



£2 = CD{^J)OD + cso{^J)o, 



so 



(2) 
(3) 



The coefficients of the operators in Eq. (j2)) are known in the leading logarithmic approxi- 
mation 1^ [n I 



Cfc(/i) 



dsirn) 



-3//3o 



Ce(/^) 



Ois{m) 



-8/(3/3o) 



(4) 



where /5o = 11 — |ri/ is the one-loop coefficient of the beta function for nf active quarks. The 
coefficient of the last operator in Q is gauge- dependent; the value quoted above corresponds 
to the Feynman gauge. 

To order 1/m^ there are two local operators contributing at tree- level, the Darwin term 
and the spin-orbit interaction energy respectively 



^9. 



(5) 
(6) 

The coefficients Ci5(/i) and csoif^) were calculated in (see also Q,|l3) 

using a different 

operator basis and in the absence of light quarks. We reconsider here their calculation in a 
basis more suited for our present purposes, with a different result for Cd{ij)- 



In addition to the operators shown in 
operators which must be added to C2: 

a) the gluonic operators 



there are three other types of dimension-6 



1 

4' 



- -9fabcF'^yFl)^F{^ , 



a\v\ 



(7) 



A third gluonic operator F'^^D'^F"'^^ has been eliminated in zero momentum insertions 
with the help of the Bianchi identity, which gives 

j d^x(40^3 + O^DF). + -^F^^D^F'^^'') = (8) 

Their coefficients have been computed many times |l^ ll^ • At one-loop order they 
are given by 

as(m) , , 2a., fm 



Ci?3(m) 



457r 



C(DF)2(m) 



ISvr 



(9) 



At one-loop order Ops does not mix with other dimension-6 operators and will be 
neglected in the following. The second operator 0(^of)^ will be eliminated in favor of 
4-quark operators using the equation of motion of the gluon field (see Eq. below). 

b) Operators which vanish by the equation of motion of the heavy quark field iv ■ Dh = 0. 
Even though their expectation values vanish, they can contribute when considering 
mass corrections to the matrix elements of currents. A discussion of their renormal- 
ization is given in Appendix A. 

c) 4-quark operators built out of heavy-light and light-light quark fields. It is these 
operators which will be the main point of interest of this paper. 

There are four independent heavy-light 4-quark operators, which can be conveniently 
chosen as follows 

^i' = yE(^"V'^^)(^7'^^"^) (10) 

Q 

^2 



9_ 
2 

Q 

^2 



(^2 = '-:^Y.^qi,ir<l){hrirh) (11) 



= ^-^Y.^^l,^)^h'h) (12) 
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2 

2 



= -:^Y.^<ll,im){hl'l^h). (13) 

The summation extends over nj dynamic quarks, for example rij = 4 (g = u, d, s, c) for the 
case of a heavy bottom quark. 

The equation of motion for the gluon field gives a relation among the Darwin term Op, 
Eq. and one of the 4-quark operators Of' 

OD = Of + ^{ht''h){hrh). (14) 

The structure of the possible 4-quark operators containing only light quark fields is similar 
to that of the heavy-light operators. There are altogether 10 such operators, which can be 
divided into three groups with different transformation properties under SU {n f) SU {n f) r. 
The operators in each group mix only among themselves under renormalization. Their 



renormalization has been discussed in detail in Refs. I2G 



21 



Only one group of operators 



4 



is relevant for our case, which can be chosen as 

-.2 



^1 = |-E(^v^^)(^7^^"?) 
<i,<i' 



(15) 
(16) 
(17) 
(18) 
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The complete basis of the dimension-6 operators includes also three nonlocal operators 
consisting of time-ordered products of dimension-5 operators 



o 



d^xT[h{iDfh]{x) [h{iDfh]{id) 
A'^xi:[h{iDfh]{x) [/i^a-F/i](0) 



J d'xT[h^-a-Fh]ix) [h^-a-Fhrn. 



(19) 
(20) 
(21) 



We will use a compact vector notation for the thirteen operators (jSEEll^)) defined as 
follows 



O 



hi 



o 



with Ti 



Oso 
Okk 



o 



hi 



Of 

hi 



\ Of J 



(22) 



\oil 



and a corresponding vector notation for their Wilson coefficients Ci{fi), carrying the same 
index as the operators Oi. Requiring that the total renormalized Lagrangian C2 be scale- 
independent gives a renormalization-group equation for the coefficients c(/i) 



d 



/i— c — 7^c = . 
dfj, 

The anomalous dimension matrix can be written as 



(23) 
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^ A B 0^ 
C D 
E F 



(24) 



To one-loop order the blocks in this matrix take the values 

/ 4 \ 



A 



9' 



(4vr) 







146 
9 





-12 6 

\ -10 12 y 



B 



f 


9 








\ 







-5 





8 
3 




-18 
















40 
3 





64 
9 


\ 


9 


5 
3 





9 / 



c 



/|n/ \ 





\ 0/ 



D 



9^ 



{An) 



, . 4 4 4 4 , 

jdiag(13 - -Uf, 13 - -Uf, 22 - -nj, 22 - -n/J 



+ \ 



E 



(4 



/ 113 



9' 



(An) 



71 



+ 1^/ 



v 



9 ' 3'*/ ^ 

f 13 - Inj 
I 12 








/ 





3 






1 44 
\ 3 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



22 - 

22 - ) 

The blocks D and F contain contributions which take into account the running of the g 
factors contained in the definition of the four-quark operators. We have implicitly used the 



fact established in 



13| that the form of Oso as given in ® is preserved under renormaliza- 



tion. We agree with [l5| on the renormalization of the Darwin term in A in the absence of 
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the mixing with the four-quark operators. The lower diagonal block F has been calculated 
previously in ISl. {Note added. After correcting the entries Cn, Dn, En we agree with the 
results of Ref. j30|. These changes affect only the result for the Darwin term C£)(/i), but 
leave all the other results of the paper unchanged.) 



III. MATCHING AND RUNNING 

The coefficients of the operators Ti. in Eq. ()22|1 are given, at the matching scale /i = 
m, by their tree-level values c-}^{m) = column(l, 1, 1, — 1, 1). The coefficients of the non- 
local operators ()19H2H) are given simply by the products of the respective local dimension-5 
operators 

Cfcfc(/i) = 1 (31) 
CmM = ( . (33) 



Integration of the RGE for the Wilson coefficient of the Darwin operator CD(/i) is compli- 
cated by the mixing between Od and O^K The initial conditions for the leading log running 
of the Wilson coefficients of the 4-quark operators are cf'(m) = 0. We compute below also 
the leading nonvanishing contribution to the matching of these Wilson coefficients, of order 
a,{m) (see Eqs. fP^ -dHHl)). 

We start by solving the RGE for the coefficients C{fj,) = (cd(/x), c5''(/x)), which has the 
form 

/x|-(7(/x) = AlCif,) + 5(/i) (34) 

with A2 the corresponding 2x2 block of the anomalous dimension matrix 7 in Eq. 
and -B(/i) an inhomogeneous term 

A--^(' 5% g,. ^^{-f- 1°--"'° I (35) 

' 4" 9 13-1,,, ' 4i _i8 + 9,-6/ft I ^ ' 
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We denoted here z = as{^j) / ois{m)-, and we neglected contributions to B{^) from the light- 
light 4-quark operators Of. These contributions are of the form ^^(as log(/i/m))", and 
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are thus formally of higher order than the leading logs computed here. They are however 
calculable, see below. 

The matrix A2 is diagonalized as 

--(-r-;)' 

with 7i = 13,72 = 4 — ^nj and R = 1 + -fjUj. Redefining the Wilson coefficients as 
C'(yu) = UC'{fi), the RGEs for the components of C"(/i) are decoupled 

l^-^C'i^i) = U-'A^UC'ii^) + U-^B{ij) (37) 

and can be solved by simple integration. 

Solving the RGE for the Wilson coefficients cc^fi), c^^fi), one has to distinguish the two 
cases i) 7^ 3, and ii) Uf = 3 separately, corresponding to 72 being nonvanishing, and 
vanishing, respectively. 

The solutions for the case i) Uf ^ 3 are given by 

If 16n) + 96n. - 81 „ ,^ 1232 ^ s , , 

R \ 108(n/ + 6) 3159 ^ ^ ^ 

12^2 + 353nf + 1605 _^^/(2^,,) 77(16n} - 48n/ - 1053 1 
^ 12(n/-3)(n/ + 6) ^ ^ 3159(n/ - 3) J 

cl^f,) = - l^!ll±^z-m. + ^,-7.7(2/30) (39) 
^ ^'^^ i?\4(n/ + 6) 117 ^ ^ 

12^2 + 353n^ + 1605 77(4n^ - 51 



-z 



Uinj - 3)inj + Q) lir(n; - 3) 

The corresponding results for case ii) Uf = 3 are 

1 f 13 ^_eM._ 1^^-13/(2,0) , 154, , 94911 
1053 

1 1 13 _fi/fl„ 308 _iq//ooN 154, 725 



^'^^ i?\36 1053 9po ^ 4212 J ^ ^ 



>^(^) = ± ) ^^-6//3o _ -13/(2/30) _ — log^ + I (41) 

R\ 12 117 9po 468 J ^ ^ 

For both cases i) and ii), the sum of the two Wilson coefficients is given by a common 
expression 

cd(^) + cf (/.) = - ^.-3/(2,0) + 308 _ ^^2) 

This combination corresponds to the coefficient of the Darwin operator in Ref. in a 
basis in which the 4-quark operator is eliminated using the equation of motion Eq. ()14|) . 
The result Eq. ^ agrees with Eq. (19) of Ref. ^^\. 



The results for cz)(/i), cW^) presented above neglect mixing into O^^ from the light-light 
4-quark operators Of . As mentioned above, they are formally next-to- leading order, and are 



ly computed. They require the 
23|), and the 2- loop anomalous 



comparable with other contributions which can not be yet fu 
one-loop matching of the Darwin operator (obtained in Ref. 
dimension (not yet available). However, the mixing of cf into CD(/i), c5*'(/i) is a well 
defined effects and can be computed from our results. 

We present next the Wilson coefficients cf in leading log approximation. The equation 
of motion for the gluon field relates the operator 0(^ofY ^'^ 4-quark operators 

O^^DFY = \ (D^Pp (^aF'^^") = 20f + Of + ^Cht^h) {heh) . (43) 

We will use this relation to eliminate the gluonic operator in favor of the 4-quark operators 
Of' and Of. Using Eq. Q, this gives the initial condition for cf{m) 

11 , ^ 2as(m) , , 

cf{m) = -^^. (44) 

It is now straightforward to solve the RGE for the coefficients cf{fi). We focus on the 
physically interesting case of the corrections, and take as thus Uf = 3. The results are 
given by (with (3q = 9) 



(mj 
(m) 



cfifi) = (0.0246z"'^^/(^*) + 0.43592~'^2/(2^») + 0.16462^^^/(2^") + 0.37A8z~^^/^^^''^)cf{ 
cf{^) = (-0.0434z"^i/(2/3o) ^ o.l96l2-'^2/(2^«) - 0.37642-^«/(2/3o) ^ o.2237z-'^*/(2'^"))c'/( 
4'(/i) = (-0.0512^^1/(2^°) ^ 0.07672"^2/(2A') + 0.074l2~'^3/(2'3o) _ 0.09982"^*/(2^"))c'/(m) 
4'(/i) = (0.02902-''i/(2'3o) + 0.17062~'^2/(2^«) - 0.0324z-^3/(2'3o) _ 0.16722"'^*/(2*))cf (m) 

(45) 

The eigenvalues of the mixing matrix F are given by a;s/(47r)(5j, with 5i = 20.2680,52 = 
24.8138,53 = 4.4516,^4 = 12.0222. 

We turn next to the spin symmetry breaking operators. We reproduce the result for the 
Wilson coefficient of the spin-orbit operator Oso obtained in Ref. [l^ 

csoi^^) = 2 ( '^'°-l. (46) 



as(m)_ 

This coefficient has been also obtained without an explicit calculation in uM from 

. No such RPI constraints exist for the 



coefficient of the Darwin term cd, contrary to the argument of Ref. Q|. A discussion of the 
RPI constraints for the renormahzed HQET Lagrangian is given in Appendix B. 

Finally, the coefficients of the remaining heavy-light four-quark operators (|11II13|) are 
found to be 



c-(/x) 



10 



.^-3//3o + 



.-6//3o 



15 



+ 



21 -Auf 3 - Auf 

20(99 + 96n/ - 16n}) 
(4n;-39)(4ny-21)(4n;-3) 



3(12 -n 



.^-3//3o + 



39 - Auf 

^-(13-|n/)/(2/3o) 
.^-6//3o 1 



+ 



f) 

4f-639 + 156n 



3 15-2n 



f) 



33 - 2n 



f 



(47) 

(48) 
(49) 



/ 



3(n/ - 12)(2r2j - 33)(2r2j - 15) 

A full one-loop determination of the coefficients cYi^Ji) requires an explicit matching 
calculation. This involves computing the box diagrams for heavy-light quark scattering in 
QCD shown in Fig.l. When expanded in powers of 1/m the total result for these diagrams 
is 



'QCD 



(4vr) 



2A 



27A 



+n (i^::;2 - 2 ) [7^75] /3a [7^75] 57 



+ 0(l/m=^ 



(50) 
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To simplify the calculation we have taken a massless light quark scattering in the forward 
direction. The infrared singularities have been regulated with a finite gluon mass A. 

A computation of the same diagrams using the effective theory Feynman rules for the 



heavy line with insertions of all 1/m corrections up to second order gives, in the MS scheme, 
,4 



I- HQET 



(47r) 



2A 



m 



-^b.f'Wb'n, + ^[7.75t1/3«[7'^75t1.. 



Stt 



27A 



[7M75]/3a[7''75]57 



+ 



^ (^n^ + 2]b,l5t''WbV]s^ 



2 / ;\2 
+ 9(^l^;;^ + 2)[7/.75]/3a[7^75]57 



+ 0(l/m3 



(51) 



The ffist two terms in the 1/m expansion agree identically, as does the logarithmic de- 
pendence on the infrared regulator A in the 1/rn? part. Imposing equality of these two 
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expressions requires adding the four-quark operators p0lll3p to the effective theory with 
coefficients 



18l_ 
(47r)2 

(47r)2 


(47r)2 



10 fi' 
^log — 

6 m 



16 

— log — 



40 

y 

64 



(52) 

(53) 
(54) 
(55) 



9 m 

The logarithmic terms agree, as they should, with the leading terms obtained by expanding 
the LLA expressions Eqs. 

One should point out that the constant terms in the full one-loop expressions for 
these coefficients are scheme-dependent. The 0{as) terms become well-defined only af- 
ter renormalization-group improvement. In this process the renormalization scheme depen- 
dence of the constant terms in ()52II55|) is cancelled by the scheme dependence of a two-loop 
anomalous dimension [2^. However, there are reasons to expect the partial results ()52l|55 1 
to provide a good estimate of the full 0{as) correction, based on empirical evidence [25 1 
that the contribution of the two-loop anomalous dimension is often much smaller than the 
constant term in the one-loop expression of the coefficient. 

On the other hand, the logarithmic terms in ()52II55|1 and the LLA sums ()39p47Ml?|l are 
scheme independent. However, the use of the leading- log resummed expression is appropriate 
only when the logarithmic terms in ()52II55|) dominate clearly over the constant terms. If this 



a 



P 



a 



P 



Y by 8 

FIG. 1: Diagrams contributing to heavy-light quark scattering which are needed for the matching 
conditions for the four-quark operators c^'(m). The thick line denotes a heavy quark and the 
dotted lines denote gluons. 
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is not the case, then the use of the full one- loop results ()52II55|) . with all the above limitations, 
is likely to be a better approximation to the true value of the Wilson coefficient. Such a 
situation is encountered in the case of the charm quark heavy mass expansion, as discussed 
in the next section. 



IV. APPLICATION: HEAVY MESONS HYPERFINE SPLITTING 

The importance of considering the four-quark operators ()10II13|) becomes immediately 
apparent if one notes that the equation of motion for the gluon field gives Od = O^K The 
coefficients of these two operators (|38|) and (|39p get contributions of the same order, such 
that neglecting O^^ would lead to an incomplete result. Adding the contribution of 
together with that of the Darwin term gives, in the leading logarithmic approximation, the 
"effective" coefficient Eq. (jl^ for the local operator of 0{l/m?) 



, . , 30^ 



as{m) J 117 



(56) 



It is possible to eliminate Oi'' from the effective Lagrangian to all orders in as, in favor of 
the Darwin operator, which contributes with coefficient C£)(/i) + cj''(/i) — > C£)(yu) [30]. 



More interesting from a practical point of view is the fact that the operators O!^^ and O^'' 



break the heavy quark spin symmetry and thus contribute to the hyperfine splittings of the 
heavy hadrons. The remainder of this Section is dedicated to a discussion of these effects. 

We first briefly review the derivation of the mass formula for the mass of a heavy hadron, 
including the contributions of the 4-quark operators (fT0|l - |T3|l . The mass of a heavy hadron 
can be written as an expansion in inverse powers of the heavy quark mass as 5|, |6[ 

uibm = mb + A+ {THint exp ^-i j d^xHintix)^ ) (57) 

with Tiint = —j^int = ~2m^'^i — ■ ■ • the interaction Hamiltonian. The expectation value 
is taken between eigenstates of the leading order HQET Lagrangian which have the mass 
rrib + A. This gives, including terms up to order l/ml, 

m^M = + A - ^^{Ci) - ((A) + 2{cMy{Okk) (58) 

+ Ckifi)cM{Okm) + 2{cMf{Omm) + (/i) (Of ) ] + • • • . 



i=l 
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The expectation values appearing in this formula can be parametrized as (using the notations 
of Ref. B) 

(Ci) = -Ckifi)f4- cUf^)^nl (59) 

(£2) = Cz5(/i)/if) + y^C5o(/i)/i|o (60) 

(Okk) = /iL (61) 
{Okm) = ^f^lo (62) 

(Omm) = f^GG + I^pI + \dHp\ , (63) 

with dn = for a pseudoscalar meson and -1 for a vector meson. 

The matrix elements of the four-quark operators (fTn|) - ()13j) can be estimated with the help 
of the factorization approximation 3]. (A n alternative model-independent determination 
of Of"' (based on a method proposed in has been given in Ref. j3].) This is done by 
first applying a Fierz transformation to bring the operator into the form {hT'q){qT'h), after 
which the vacuum state is inserted between the two currents. Only color singlet currents 
will give a nonvanishing contribution. In this way we obtain 

2^^(5(1^)1^ 5^(g7.75t'^9)(/^7^5t'^/^)|5(t;)) = (64) 

,2 



^(5(t;)||-5^(g7.75g)(/^7'^75/^)|S(t;)) = -^/>B (65) 



q 

,2 



^(5*(t',5)||-5^(g7.75t'^g)(/^7'^75t'^/^)|5*(t;,e)) = ^/Im^ (66) 



9 

^2 



^(5*(t;,e)||-5^(g7.75g)(/^7^5/^)|5*(^;,5)) = ^/Im^ • (67) 

We used here the following relations which are valid in the heavy mass limit 

{B{v)\h'yf,'y5q\0) = ifsmEVf, (68) 

{B*{v,e)\ht.q\0) = fBmBe;. (69) 

Combining these relations one obtains the following expression for the hyperfine splitting 

rriB* - TUB = + 3(2^^)2 {(^so{p)plo + cMcMpIg (70) 
HcMfpl - 4\p)^fBmB - c'Ap)'^flmB) + 0{l/ml) . 
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LLA 


One-loop 


II = 0.5 GeV 


= 1 GeV 


= 0.5 GeV 


II = 1 GeV 


c'i\iJL/mc) 


-0.268 


-0.074 


-0.653 


-0.503 


cfil^/rric) 


-0.117 


-0.038 


-0.348 


-0.268 




-0.401 


-0.232 


-0.577 


-0.483 




-0.156 


-0.105 


-0.308 


-0.268 



Table 1. The coefficients of the spin-symmetry violating four-quark operators 
for the charm and bottom cases, at two values of the factorization scale n — 0.5 
GeV and = 1 GeV. 

This relation can be written only in terms of observable quantities by expressing the 
quark mass in terms of the pseudoscalar meson mass and multiplying with uib + niB* 

m|* - m| = -CmilA"^ + ^(^2mB) + cso{fJ')fJ'SO + cMcmifJ^fJ'la (71) 

+{cm{n)rp'A - 4\^^f-^fBmB - cWiif-^flruB^ + 0(l/m|) . 

In Table 1 we present values for the coefficients of the two four-quark operators which 
contribute to the hyperfine splittings. They are given in both the leading-log approximation 
and the "full" one-loop form at two different values of the factorization scale // = 0.5 
and 1 GeV. In computing these values we used h-qcD — 250 MeV and rric — 1.39 GeV, 
mb — 4.8 GeV. For the bottom hadron case we neglected the change in the running of ag 
across the charm threshold, which gives a negligible error. The one-loop result suggests that 
the leading-log approximation is very poor for the charm case, where the logarithmic term 
accounts for at most 30% of the total correction. The situation is slightly better in the 
bottom quarks, where its contribution is enhanced to about 50%. 

For the purpose of illustration we show in Table 2 the combined contributions of these 
two operators to the hyperfine splitting of the D and B mesons using for their coefficients 
both approximation methods. For the reasons discussed above we tend to favor the full 
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one-loop results over those obtained from the leading log approximation. In computing the 
matrix elements of the operators in Table 2 we have used the hadronic parameters = 170 
MeV, fe = 180 MeV, = 1.97 GeV, = 5.28 GeV. 







W)d'-d 




{Of)B*-B 


Afl (MeV) 


Ab (MeV) 


0.5 GeV 


(543 MeV) 3 


(493 MeV)3 


(783 MeV)3 


(712 MeV)3 


7.4(LLA) 
18.9(1-1.) 


2.7(LLA) 
4.2(1-1.) 


1.0 GeV 


(431 MeV)3 


(392 MeVf 


(622 MeVf 


(565 MeV)3 


l.l(LLA) 
7.3(1-1.) 


0.8(LLA) 
1.8(1-1.) 



Table 2. Matrix elements of the four-quark operators O2 and in the factor- 
ization approximation for two different values of the factorization scale /i = 0.5 
and 1 GeV. In the last two columns the total contribution of these two operators 
to the hyperfine splittings of the D and B mesons Aj = tuh* — is shown in 
both the leading-log approximation and using the one-loop result (1-1.) for their 
coefficients. 

The relation ()71|) contains the as yet unknown parameters f^^soy^^nG P\- '^^^ ^^^^ 
parameter vanishes exactly for S-wave mesons in potential models as the corresponding 
operator © corresponds to the spin-orbit interaction energy of the heavy quark p]. The 
remaining two parameters are not easy to estimate in a model-independent way. The con- 
tributions of the four-quark operators are shown in the last two columns of Table 2. Their 
total contribution is positive and amounts to about 10-20 MeV in the charm case and 1-5 
MeV in the bottom case. 

Comparing with the measured hyperfine splittings one can see that a positive l/m^ 
correction to (f7T|) agrees with the sign of the correction to the heavy mass scaling law. 
Unfortunately, the lack of information on the precise values of /il^, /U^^., p\ prevents us from 
making a more quantitative analysis of the deviations from leading order scaling. 

Finally we note the relevance of these corrections for the model-independent determi- 
nation of the matrix element /i^ from the measured hyperfine splittings in the B system. 
From (f7T| we obtain 0.358 GeV^ = /Xg -|- 0.026 -|- 0.021 + terms proportional to fJ^so/^B 
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and fi^Q/niB- The two numbers on the right-hand side correspond to A/i^ and respectively 
to the four-quark operators' contribution. In this estimate we used A = 400 MeV and fiQ 
=0.35 GeV^. This shows that the theoretical error of this determination coming from the 
l/ms terms in (j71l) may be as high as 13%. 
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APPENDIX A 



A complete study of the renormalized HQET Lagrangian at order requires also the 
computation of the coefficients of operators which vanish by the equation of motion of the 
heavy quark field. There are altogether six such operators: three local operators which will 
be chosen as in j3| 

Os = ^ha^,F^'{v ■ D)h + ■ D)a^,F^'h (72) 

O9 = ihD^{v D)h + ih{v- D)D^h (73) 
Oio = hiivDfh, (74) 

and three nonlocal operators consisting of time-ordered products of dimension 5 operators 

O^^ = ^-j d'^xTlhiiv ■ Dfh]{x)[h{iv Dfh]{{)) (75) 
Oke = i j d'^xT[h{iD fh]{x)[hiivDfh]iO) (76) 
Ome = i j d^xT[h^a ■ Fh]{x)[h{iv Dfh]{Q) . (77) 

They will be combined together in a vector Tig defined as 

He = C0lumn(08 , O9 , ^lO , O^e , Oke , Ome) • (78) 

These operators mix only with the operators in Ti. fl22|l . The renormalized operators Ti 
and Tie satisfy a renormalization group equation which can be written in matrix form as 



. (79) 
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The block A has been given in ()25p. and the remaining ones are given by 
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(80) 



(81) 



At the scale /i = m the coefficients of the operators in (|75|) are given by 

Ce = column(0 ,0,0,1,-1,1). 



(82) 



Integrating the RG equation (f?^ with this initial condition one obtains the following expres- 
sions for the coefficients of the three local operators in (|78p in the leading log approximation 

9 



4 4 2po 

C9(/i) = 3;2-s/(3/3o) _ 3 
3 



cio(/i) 



-6//3o _ ^ -2//3o , g;2~16/(3/3o) , ^^-8/(3/3o) _ ^ 

9 30 18 



(83) 
(84) 
(85) 



We denoted here z = {ag^jj) / oisirn)). In deriving (|H^ we used the result derived in Sec. Ill 
for the Wilson coefficient of the Darwin term. It is important to emphasize that these coeffi- 
cients are (just as Ce(/i) in Q) gauge dependent. Their expressions given above correspond 
to the Feynman gauge. 
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APPENDIX B 



The matchin g re lation connecting the heavy quark field Q in QCD with the effective 



theory field h is 



3 



Q{x) = Qexp{—imv ■ X'^)h{x) 



(86) 



with 



Q = exp 



exp 



(87) 



2m J ^ yAm? 

and p\\ = v-D-^, p± = p — p\\. The field h contains both "upper" and "lower" components 
h±, satisfying ■f>h± = ±h±. One has 



exp(— -imf ■ x'^)h 



-h + e' 



1 



-h 



(88) 



It can be shown jl^ that the field transformation (jHTjl decouples the two components /i+ 
and h-. The HQET Lagrangian (PE)) used in the main text refers only to the /14. part. 

The most general form for the 1 / term in the HQET Lagrangian includes, in addition 
to the operators introduced so far, also the operators 



'-lha,,F^'''{v ■ D)h + |/i(i; ■ D)a^,F^^h 



Og = ihD^{v- D)h + ih{v D)D^h 
do = hiivDfh. 



(89) 
(90) 
(91) 



Their renormalization properties are studied in the Appendix A. 

The purpose of this Appendix is to show that requiring the invariance of the HQET 
Lagrangian ^ under a small change in the velocity v (the so-called reparametrization in- 

n. 



variance [16*) fixes the coefficients csoiu. 



is obtained in this way, in contrast to 



[14i] and Cg{fi). However, no constraint for Cd{ij) 



We start by computing the change in the /i+ field under an infinitesimal change of the 
velocity parameter v ^ v' = v + Av. This can be obtained from (j86|) by taking into account 
the invariance of the QCD field under this transformation 6Q = as 

'1 + ^ 



5h, 



S 



■ exp{imv ■ x)Q ^ 



exp{—imv ■ xf)h . 



(92) 
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Explicitly, to first order in Af and up to second order in 1/m, 



(93) 



Note that no negative component field h- is introduced by this transformation. 

The variation of the HQET Lagrangian (POj) is obtained by the variation of the effective 
field dnSj) and the variation of the f-dependent operators, including the operators 
After some algebra one obtains 



5C = (1 - Ck{^i))h+iAv ■ Dh+ - — {-(ce + cg - l)h+{Av ■ D){v ■ D)h^ 

cso + 1~ 



^ I Ce + 1 + Cg + Cm + 2 



- ^Ce + 1 + Cg - Cm - j K^gYl ■ Avv>^F,,h+0{l/m^] 

This will vanish for any Av provided that the following identities hold 



(94) 



Cfc(/i) = 1 

= -1 - Ce(/i) 
C5o(/^) = -1 - 2Cm(yU) . 



(95) 
(96) 
(97) 



The first constraint has been given in 16| and is perhaps the best known application of 



the reparametrization invariance principle. The constraint ()97p has been presented in 
and its prediction for cso agrees with the exphcit calculation in [l^ (see (jlBj) ). The relation 
(jnni) together with (jH) predicts the following value for cg(/i) 



C9(/i) = 3 



Q^(/i) 



-8/(3A)) 



1 



(98) 



This agrees with the expression of this coefficient obtained by direct computation in leading 
log approximation 
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